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The effects of different forms of the sound- velocity function c s (T) on the hydrodynamic evolution 
of matter created in the central region of ultra-relativistic heavy-ion collisions are studied. At high 
temperatures (above the critical temperature T c ) we use the sound velocity function obtained from 
the recent lattice simulations of QCD, whereas at low temperatures we use the ideal hadron gas 
model. At moderate temperatures different interpolations between those two results are employed. 
They are characterized by different values of the local maximum (at T = 0.4 T c ) and local minimum 
(at T = T c ). The extreme values are chosen in such a way that at high temperature all considered 
sound-velocity functions yield the entropy density consistent with the lattice simulations of QCD. 
We find that the presence of a distinct minimum of the sound velocity leads to a very long (~ 20 
fm/c) evolution time of the system. Since such long evolution times are not compatible with the 
recent estimates based on the HBT interferometry, we conclude that the hydrodynamic description 
becomes adequate if the QCD cross-over phase transition renders the smooth temperature variations 
of the sound velocity, with a possible shallow minimum at T c where the values of c 2 B (T) remain well 
above 0.1. 
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I. INTRODUCTION 

The analysis of the data collected recently at RHIC 
suggests that matter created in ultra-relativistic heavy- 
ion collisions behaves like a perfect fluid [l|. This hy- 
pothesis inspires new interests in the studies of rel- 
ativistic hydrodynamics of perfect and viscous fluids 

S, H H 0, 0, S 0, OH EI Ei Ei EI Ei EE E3 . The 

aim of this paper is to study the effects of the temperature 
dependence of the sound velocity on the hydrodynamic 
evolution of matter created in ultra-relativistic heavy-ion 
collisions. We use the approach to the relativistic hydro- 
dynamics of perfect fluid where the sound velocity is the 
only thermodynamic parameter entering the formalism 
[IsL EH 01 • In this way, we are able to observe and ana- 
lyze the immediate consequences of a particular choice of 
the function c s (T) on the space-time evolution of matter. 

Our study may be regarded as complementary to Ref. 
[2l| . where the effects of a particular choice of the equa- 
tion of state on the particle spectra and elliptic flow were 
studied. Our approach differs from Ref. [2l| by the se- 
lection of the equations of state under investigations. In 
our case the equations of state are always constrained 
by the lattice results in the high temperature limit and 
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differ from each other mainly in the neighborhood of the 
critical temperature. On the other hand, the main differ- 
ences between the equations of state studied in Ref. [2l| 
appear in the region above the critical temperature. 

In our approach the full information about the equa- 
tion of state is contained in the temperature dependent 
sound velocity. At high temperatures, T > 1.15 T c , we 
use the sound velocity function obtained from the recent 
lattice simulations of QCD (22y, whereas at low tem- 
peratures, T < 0.15 T c , we use the result of the ideal 
hadron gas model [19]. At moderate temperatures dif- 
ferent interpolations between those two results are em- 
ployed. The interpolating functions have a local maxi- 
mum at T = 0.4 T c (corresponding to the maximal value 
of the sound velocity in the hadron gas) and a local min- 
imum at T — T c (corresponding to the expected mini- 
mal value of the sound velocity at the phase transition) . 
The values of the sound velocity at the maximum and 
minimum are chosen in such a way that in the limit 
of very high temperatures all considered sound-velocity 
functions yield the entropy density consistent with the 
lattice simulations of QCD. 

We emphasize that in all the considered cases we 
take into account only smoothly varying functions c s (T), 
hence we deal with the cross-over phase transitions rather 
than with the rigorous phase transitions of a given or- 
der that might be determined by the discontinuity of the 
appropriate thermodynamic variable. Nevertheless, the 
deeper is the minimum of the sound velocity at the crit- 
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FIG. 1: Temperature dependence of the square of the sound 
velocity at zero baryon density. The plot shows the result of 
the lattice simulations of QCD [22] (solid line) and the result 
obtained in the ideal hadron gas model [H (dashed line). A 
piece of the thick solid line describes the simplest interpolation 
between the two calculations. The critical temperature T c 
equals 170 MeV. 



ical temperature, the more the system evolution resem- 
bles the behavior typical for the exact first order phase 
transition. The main difference between the first order 
phase transition studied commonly in the literature and 
the approximate first-order phase transition studied be- 
low (as one of the examples) is that in our case the phase 
transition does not lead to the ideal quark-gluon plasma 
but rather to strongly interacting plasma characterized 
by the lattice results. 

In spite of such differences, in agreement with earlier 
studies (23l . l2~il . I25L we find that a deep minimum of 
the sound velocity leads to the very long, exceeding 20 
fm (c = 1), evolution time of the system. We further 
find that even moderate minima lead to quite long evo- 
lution times of the system, which exceed 15 fm. Clearly, 
such long evolution times are not compatible with the 
recent HBT estimates of the lifetime of the system. For 
example, one finds times of about 10 fm using th e RHIC 
data in the relation i?i(mr) = T^Tk/niT [27j], which 
connects the longitudinal pion correlation radius Rl, the 
kinetic freeze-out temperature T^, the evolution time r, 
and the transverse mass of the pion pair my. Another 
example is Ref. [28] that obtained only 6 fm for the effec- 
tive duration of the hydrodynamic evolution in Au+Au 
collisions at RHIC. 

Only if the sound velocity has a shallow minimum re- 
sulting from the simplest interpolation between the ideal 
hadron gas result at T = 0.85 T c and the lattice result at 
T = 1.15 T c , see Fig. [TJ one finds that the evolution time 
of the system might be compatible with the HBT result. 
We thus conclude that the hydrodynamic description be- 
comes adequate if the QCD cross-over phase transition 
renders the smooth temperature variations of the sound 




FIG. 2: Three different forms of the sound velocity used in the 
present hydrodynamic calculations. The solid line describes 
the interpolation between the lattice and the hadron-gas re- 
sults [l^ ] with a shallow minimum where c 2 s — 0.14 (case I), 
the long-dashed line describes the interpolation with a dip 
where ci = 0.08 (case II), finally the dashed line describes 
the interpolation with a deep minimum where (? s = 0.03 (case 
III). In order to render the same entropy density at high tem- 
peratures, in the cases II and III the sound velocity in the 
region T « 0.4 T c is slightly increased compared to that con- 
sidered in the case I. The larger values of the sound velocity in 
this region may be attributed to the repulsive van der Waals 
effects. 



velocity, with a possible shallow minimum at T c where 
the values of c^(T) remain well above 0.1. 



II. TEMPERATURE DEPENDENT SOUND 
VELOCITY 

In Fig. [TJwe show the function c 2 s (T) obtained from the 
ideal hadron gas calculation and the lattice simula- 
tions of QCD [12]. The ideal hadron gas result (dashed 
line) was obtained with the full mass spectrum of reso- 
nances [1^, [3(3] ■ One can check that at very small tem- 
peratures the sound velocity is determined basically by 
the properties of the massive pion gas; for T « m w , 
where is the pion mass, one finds c s (T) = ^jT/m^. 
Clearly, this behavior is quite different from the limit- 
ing value c s = 1/V3 characterizing the massless pion 
gas. The lattice calculations (solid line) were obtained 
for physical masses of the light quarks and the strange 
quark [221 - O ne can see that the results of the two calcu- 
lations cross around the critical temperature, we assume 
T c = 170 MeV, and may be naturally interpolated (see a 
thick solid line in this region) . Taking the lattice result at 
the face value, one expects that the sound velocity signif- 
icantly drops down in the region T » T c . Similar behav- 
ior, with c s (T c ) reaching zero, is expected in the case of 
the first order first transition where the changes of the en- 
ergy happen at constant pressure. However, the lattice 
simulations suggest that for three massive quarks with 
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FIG. 3: The temperature dependence of the entropy density and energy density, panels a) and b), as well as the energy density 
dependence of the pressure and sound velocity, panels c) and d). One can observe that the deeper is the minimum of the sound 
velocity function, the steeper is the increase of the entropy density and the energy density. 



realistic masses we deal with the cross-over rather than 
with the first order phase transition, hence the sound ve- 
locity remains finite, as is consistently shown in Fig. [TJ 
Nevertheless, the exact values of the sound velocity in 
the region T ~ T c are in our opinion poorly known, since 
the lattice calculations are not very much reliable for 
T <T C and, at the same time, the use of the hadron gas 
model with vacuum parameters becomes unrealistic for 
large densities (temperatures) [3{J. In this situation, it is 
practical to consider different interpolations between the 
lattice and hadron-gas results and to analyze the physical 
effects of a particular choice of the interpolating function. 

In this paper we consider three different sound- velocity 
functions c s (T). Below, we refer to these three options as 
to the cases I, II and III, see Fig. [21 In the case I, we use 
the sound-velocity function which agrees with the ideal 
hadron gas model of Ref. [l^| in the temperature range 
< T < 0.85 T c and with the lattice result in the temper- 
ature range T > 1.15 T c . In the region close to the critical 
temperature, 0.85 T c < T < 1.15 T c , a simple interpola- 
tion between the two results is used. We have checked 



that such a simple interpolation yields directly the en- 
tropy density consistent with the lattice result. Namely, 
the use of the thermodynamic relation 



s(T) = s(T )exp 



dT 



To 



(1) 



relating the entropy density with the sound velocity for 
zero baryon chemical potential, gives the function s(T) 
which agrees with the lattice result at high temperatures, 
s(T)/T 3 » 12 at T= 1.5 T c (H- 

In the cases II and III, the sound- velocity interpolating 
functions have a distinct minimum at T — T c . Compar- 
ing to the case I [with c s {T c ) = 0.37 and c 2 s (T c ) = 0.14], 
the value of the sound velocity at T = T c is re- 
duced by 25 % in the case II [where c s (T c ) = 0.28 
and c 2 s {T c ) — 0.08], and by 50% in the case III [where 
c s (T c ) = 0.19 and c 2 s (T c ) = 0.03]. From Eq. © one con- 
cludes that the relative decrease of the sound velocity at 
T c leads to the relative increase of entropy density for 
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high temperatures. Hence, in order to have the same 
value of the entropy density at high temperatures, a de- 
crease of the sound velocity function in the region T~T C 
should be compensated by its increase in a different tem- 
perature range. For our interpolating functions in the 
cases II and III we assume that the values of c s (T) in 
the range 0.15 T c < T < 0.85 T c are slightly higher than 
in the case I, see Fig. [2l Such modifications may be re- 
garded as the parameterization of the repulsive van der 
Waals forces in the hadron gas. The values of the max- 
ima are chosen in such a way that the entropy densities 
for three considered cases are consistent with the lattice 
result, see the upper left panel of Fig. [3] where the func- 
tions s(T)/T 3 are shown. 

We stress that in the three considered cases the values 
of c s (T) in the temperature range T c < T < 1.25 T c re- 
main significantly below the massless limit 1/ys. Such a 
limiting value is implicitly used in many hydrodynamic 
codes assuming the equation of state of non-interacting 
massless quarks and gluons for T > T c , see for example 
the extended 3+1 hydrodynamic model of Ref. [12]. We 
expect that the cooling of the central part of our system 
will be significantly slower than the cooling of the sys- 
tems containing an ideal quark-gluon plasma in the core 
and this effect reflects the non-perturbative phenomena 
taking place above T c that may be attributed to the for- 
mation of strongly interacting quark-gluon plasma. We 
also note that our study differs from the 3+1 model of 
Ref. @ where the interpolation between the equation 
of state of the ideal quark-gluon plasma above T c (not 
of the lattice QCD) and the resonance gas below T c was 
introduced. We have checked that the corrections intro- 
duced in this way differ from the non-perturbative effects 
found in the lattice simulations. Moreover, Fig. 3 of Ref. 
0| (lower left panel) indicates that the sound velocity at 
low temperatures does not drop to zero as in our case. 



III. HYDRODYNAMIC EQUATIONS 

In this Section we recapitulate the main features of 
approach to the relativistic hydrodynamics, which is the 
generalization of the method introduced by Baym et al. 
[3| and has been developed in Ref. (introduction 
of the temperature dependent sound velocity and initial 
transverse flow) and Ref. [2(J (inclusion of the cylindri- 
cal asymmetry) . We rewrite the equations in such a way 
that the sound velocity is the only thermodynamic pa- 
rameter characterizing the matter. Such reformulation of 
the hydrodynamic equations is always possible in the case 
of zero net baryon density. Since we consider the evolu- 
tion of matter formed in the central region at the largest 
RHIC energies, the approximation of zero net baryon 
density is reasonable. Thermal analysis of the ratios of 
hadron multiplicities indicates that the baryon chemical 
potential at RHIC top energies is about 30 MeV, i.e., it 
is much smaller than the corresponding temperature of 
about 170 MeV [3ll. l32|. We also restrict our considera- 



tions to the boost-invariant and cylindrically symmetric 
expansion. The assumption of boost-invariance is again 
good for the central region, while the effects of azimuthal 
asymmetry are typically small, of the order of 10%, and 
have no effects on our conclusions. 

In the case of vanishing baryon chemical potential the 
hydrodynamic equations may be written in the following 
form 



u fl d l ,(Tu v ) = d v T, 
9„ (W) = 0, 



(2) 
(3) 



where T is the temperature, s is the entropy density, and 
it' 1 = 7(l,v) is the hydrodynamic four- velocity (with 
7 = 1/Vl — v 2 ). Since temperature is the only indepen- 
dent thermodynamic variable, all other thermodynamic 
quantities can be obtained from pressure given as a func- 
tion of temperature, P = P{T). Such a function plays 
a role of the equation of state. With the help of the 
thermodynamic relations 

de = Tds, dP = sdT, e + P = Ts, (4) 

other thermodynamic quantities may be obtained. In 
addition, the equation of state allows us to calculate the 
sound velocity 



dP 



s dT 
T~ds' 



(5) 



From Eq. |(5|) one immediately finds Eq. ([TJ) discussed in 
Sect. II. The complete set of the thermodynamic quan- 
tities for the cases I, II and III is shown in Fig. [3l 

For the boost-invariant systems with cylindrical sym- 
metry Eq. J3]) and the spatial components of Eq. ([2]) 
may be rewritten as 



<91nT dhxT 



91ns 

at 



r at 

(91ns 
dr 



dr 
da 

"dt 



da da 
dt 7 dr 

da 1 v r 

or t r 



0, (6) 
0, (7) 



where a is the transverse rapidity of the fluid defined 
by the condition v r = tanha. The longitudinal compo- 
nent has the well known boost-invariant form v z = z/t 
[33T |. By introducing the potential $ (T) defined by the 
differentials 



d$ = = c s ams, 



(«) 



and by the use of the two functions A± defined by the 
formula 



A± = $ ± a, 



(9) 



Eqs. (J6j) and become 
dA ± (t,r) , 



v r ± c s dA± (t, r) 



dt 



1 ± V r C s 



1 ± V r C s 
V r l\ 

r t I 



dr 



(10) 
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We note that Eq. I|10p agrees with the formalism dis- 
cussed in if one defines a± = exp(A±). If the func- 
tions A± are known, the potential $ may be calculated 
from the formula 



$> = \ {A+ + A-), 
and the velocity is obtained from the equation 



tanh 



(A, 



A. 



(11) 



(12) 



The knowledge of the function c s {T) allows us, by the 
integration of Eq. {§J, to determine $ as a function of 
the temperature; this function will be called later <E>t- 
However, to get a closed system of equations for A + and 
we have to invert this relation and obtain T as a 
function of $; this function will be called later T$. In 
this way, the sound velocity may be expressed in terms 
of the functions A + and v4_ , 



m> ( \ (A+ + A_) 



(13) 



and Eqs. 1(10)1 may be solved numerically using the stan- 
dard methods. The only restriction of our formalism 
is the condition against the formation of shock waves 
ii 



2 , rr.dcs d 

1 — cz + c s l = — 

* * dT dT 



SCs 



> 0. 



(14) 



We have checked that this condition is fulfilled by our 
sound- velocity profiles I, II and III. 

For boost-invariant and cylindrically symmetric sys- 
tems the entropy conservation law Eq. J3]) implies that 
the following quantity is conserved in time 



2ir / dr r 1 7(t, r) s(t, r) — const. 



(15) 



For one-dimensional expansion this condition is reduced 
to the famous Bjorken relation s(t) = s(to)to/t. We have 
checked that the condition (flB")) is fulfilled in our calcula- 
tions with very high accuracy, hence, no artificial entropy 
production is present in our algorithm. 



IV. INITIAL CONDITIONS 

For symmetry reasons, both the velocity field v r and 
the temperature gradient dT/dr should vanish at r = 0. 
This condition is achieved if the functions A + and A- are 
initially determined by a single function A(r) according 
to the prescription [l|3] 

A + (t = t 0) r) = A(r), A_(t = t 0) r) = A{-r). (16) 

Our main physical assumption about the initial state is 
that the initial temperature profile is connected with the 



nucleon- nucleus thickness function Ta (t) by the following 
equation 



T(t ,r) =T s [s(t ,r)]=T s 



so 



T A {r) 



Ta(0) 



(17) 



Here T g (s) is the inverse function to the entropy density 
function s(T) and the parameter sq is the initial entropy 
at the center of the system. The idea to use Eq. JT7J) 
follows from the assumption that the initially produced 
entropy density s(to,r) is proportional to the density of 
wounded nucleons at a distance r from the collision cen- 
ter [35[- We use the value so = 70.5fm~ 3 which yields 
T(t = to, r = 0) = 2T C . We note that the functions s(T) 
and T s (s) (evaluated for the cases I, II, and III) agree if 
the temperature or entropy is sufficiently large. Hence, in 
the three considered cases the initial temperature profiles 
are practically the same in the center of the system. The 
small differences appear however if we consider larger val- 
ues of r where the temperature and entropy drops down. 
We also note the global entropy defined by Eq. (fl"5)l is 
exactly the same in the cases I, II and III. 

We recall that the thickness function is defined by the 
equation 



T A (r) =2 I dzp[yjr 2 + z 2 ) 



(18) 



where the function p(r) is the nuclear density pro- 
file given by the Woods-Saxon function with a con- 
ventional choice of the parameters: po = 0.17fm~ , 
r = (1.12A 1 / 3 - 0.86A- 1 / 3 ) fm, a = 0.54 fm, A = 197. 
We note that the initial condition lfl7|) may be included 
in the initial form of the function A(r) 



A(t = t c 



<f>T<T t 



so 



Ta{t) 
Ta(0) 



(19) 



V. RESULTS 

In Fig. |4] we show the isotherms describing the hydro- 
dynamic evolution of the system with the sound velocity 
I. The numbers at the isotherms give the values of the 
temperature in the units of the critical temperature T c . 
We observe that the center of the system cools down to 
0.8 T c after the evolution time of about 15 fm. We note 
that the hydrodynamic description should be replaced 
(around T ~ 0.8 T c ) by the model describing hadronic 
rescattering (see, e.g., Ref. [§]) whose presence addition- 
ally increases the lifetime of the system. Another option 
is to assume that freeze-out happens at high tempera- 
ture (see Refs. [IB, H3| where many physical observables 
were successfully reproduced under the assumption of a 
universal freeze-out taking place at the temperature of 
165 MeV). In the latter case the lifetime of the system 
may be identified simply with the time when the system 
passes the phase transition. In the discussed case this 
time is of about 10 fm. 
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In Fig. [5] we show the isotherms describing the hydro- 
dynamic evolution of the system with the sound velocity 
II. The initial entropy is exactly the same as in the case 
I. One can notice that the central temperature does not 
drop below 0.8 T c before 20 fm. Clearly, the dip in the 
sound velocity causes a dramatic increase of the lifetime 
of the system. 

The most striking situation is presented in Fig. [U Here 
the sound velocity with the deepest minimum is consid- 
ered. Again in this case the initial entropy of the system 
is exactly the same as in the cases I and II. One can no- 
tice that the system does not pass the phase transition 
before the considered evolution time of 20 fm. Contrary, 
even at that time the system has the tendency to expand 
more, the effect indicated by the shapes of the isotherms. 
The evolution time in the case III becomes longer than in 
the case II and of course much longer than in the case I. 
The situation depicted in Fig. [6]resembles very much the 
case of the first order phase transition, where the large 
latent heat is used to increase the volume of the system at 
constant temperature, see Fig. 3 panel b) where the step 
in the energy density as the function of the temperature 
is seen. Only when the latent heat is totally consumed, 
the center of the system starts further cooling. 

Interestingly, the case III describes how our approach 
with the temperature dependent sound velocity is capa- 
ble of imitating the phenomena present at the real phase 
transitions. We have to remember, however, that the 
approximate first-order phase transition considered here 
differs from the real first-order phase transition between 
the hadron gas and the ideal quark-gluon plasma. 
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FIG. 4: Isotherms describing the hydrodynamic evolution for 
the case I; sound velocity with a shallow minimum where 
c s (T c ) = 0.37. The numbers denote the values of the tem- 
perature in the units of the critical temperature T c = 170 
MeV. 
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FIG. 5: Same as Fig. [4] but for the case II; sound velocity 
with a moderate minimum, c s (T c ) = 0.28. 
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FIG. 6: Same as Fig. [4] but for the case III; sound velocity 
with a deep minimum, c s (T c ) = 0.19. 



VI. CONCLUSIONS 

The observation of the effects related with the possible 
presence of the softest point of the QCD equation of state 
triggered the ideas about the long- living states formed in 
ultra-relativistic heavy- ion collisions (2a. WA. I25L I26L l38| . 
In view of the RHIC data indicating a short lifetime of the 
system we may analyze the problem of the lifetime in the 
reverse order, asking the question how soft the equation 
of the state may be to allow for a hydrodynamic descrip- 
tion consistent with the HBT results. In our approach the 
issue of the softest point is translated into the problem 
of realistic behavior of the sound velocity in the vicinity 
of the critical temperature. We argue that the smooth 
behavior of the sound velocity without any distinct min- 
imum is favored if we demand the short hydrodynamic 
evolution time of about 10 fm. Our calculations show, as 
expected, that the hydrodynamic evolution is very sen- 
sitive to the actual values of the sound velocity and in- 
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dicates how extremely important is the good knowledge 
of the function c s (T) . The reliable values of c s are not 
only required for the vicinity of the phase transition but 
also for smaller and larger temperatures, since the sound 
velocity plays a role of the coupling between the temper- 
ature gradients and the acceleration of the fluid. Finally, 
we note that the cooling of the system may be faster if 



the initial pre-equilibrium transverse flow is formed. The 
inclusion of such a flow in the initial conditions for the 
hydrodynamic equations leads naturally to faster expan- 
sion of matter. 

Acknowledgements: We thank Wojciech Broniowski 
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